Long range intermolecular forces in triatomic systems: 
connecting the atom-diatom and atom-atom-atom representations 
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The long-range forces that act between three atoms are analysed in both atom-diatom and atom- 
atom-atom representations. Expressions for atom-diatom dispersion coefficients are obtained in 
terms of 3-body nonadditive coefficients. The anisotropy of atom-diatom Ce dispersion coefficients 
arises primarily from nonadditive triple-dipole and quadruple-dipole forces, while pairwise-additive 
forces and nonadditive triple-dipole and dipole-dipole-quadrupole forces contribute significantly to 
atom-diatom Cg coefficients. The resulting expressions are applied to dispersion coefficients for Li + 
Li2 (triplet) and recommendations are made for the best way to obtain global triatomic potentials 
that dissociate correctly both to three separated atoms and to an atom and a diatomic molecule. 

PACS numbers: 34.20.-b,34.20.Gj,34.20.Cf,34.30.+h 
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INTRODUCTION 



There is great current interest in forming diatomic 
molecules in cold gases of alkali metal atoms, by photoas- 
sociation spectroscopy 0, 0, IE 0, by magnetic tunin; 
through atom-atom Feshbach resonances IE 0- 0- 
and through 3-body recombination [ToL In most 

current experiments, the diatomic molecules are formed 
in vibrational states near dissociation and spend most 
of their time at large internuclear separations. Once 
formed, the fate of the molecules depends on collisional 
processes which in turn depend on atom-atom and atom- 
diatom potential energy surfaces EE EE EE EE EE E3 • 

The potential energy surfaces for alkali metal trimers 
are strongly nonadditive, even for spin-polarized atoms 
(quartet electronic states) ^E^3- ^ n a ^ cases both linear 
and equilateral configurations of the M3 collision complex 
lie below the asymptotic atom-diatom energy, so that 
barrierless atom exchange reactions can occur. We have 
carried out quantum dynamics calculations including re- 
active channels for spin- pola rized Na + Na2 E3 1 Li + Li2 
0EEE3 and K + K 2 [T^. In Na + Naa, nonadditivity 
increases the well depth by about 60% yjj and increases 
the low-energy cross sections for vibrational quenching by 
a factor of 10 E2- F° r the other alkali metals the nonad- 
ditive contributions to the potential are even larger |l9| . 

Low-energy collisions are particularly sensitive to long- 
range potentials. In developing triatomic potential en- 
ergy surfaces for collision calculations, it is thus highly 
desirable to use global functional forms that include non- 
additivity and have the correct physical behaviour both 
for three-body dissociation (to separated atoms) and 
for two-body dissociation (to an atom and a diatomic 
molecule). However, although there has been a consider- 
able amount of work on both these limits, the connection 
between them has not been fully explored. This paper 
sets out to establish the connection and to suggest func- 
tional forms that are correct in both limits. 



ATOM-DIATOM LIMIT 

An atom-diatom system is conveniently described in 
terms of Jacobi coordinates R, r and 9. In this case, "long 
range" is usually interpreted to mean that the atom- 
diatom distance R is large but that the diatom internal 
distance r can be either large or small. For an atom inter- 
acting with a homonuclear diatom the long-range forces 
for R 3> r (neglecting retardation and damping) can be 
written to order R~ w j^, 

V(R,r,6) = -C 6 (r,9)R- 6 ~C 8 (r,9)R- 8 -C 10 (r,9)R- w , 

(1) 

where 



C 6 (r,9) 
C 8 (r,0) 



= C7 6 ° (r)+Cf(r)P 2 (cos 0) 
= Cg(r) + Cg(r)iMcos0) 



■C£{t)Pa (cos 9) 



= CUr)+C 2 w (r)P 2 (cos6) 
+ <# o (r)P 4 (cosd) + Cj o (r)P 6 (cos0) 



(2) 



The r-dependence of the dispersion coefficients (r) is 
central to the present work. 

Rerat and Bussery-Honvault 21] have calculated 
isotropic and anisotropic dispersion coefficients Cg(r) 
and Cg(r) for Li and Na interacting with triplet Li 2 
and Na 2 . They fitted the results to empirical func- 
tional forms based on the known asymptotic behaviour 
of diatomic polarizabilities 22]. However, they stated 
that the proper asymptotic functional form for atom- 
diatom dispersion coefficients was unknown. Merawa et 
al. [13] extended parts of this work (not including the 
r-dependence of the dispersion coefficients) to systems 
containing K and triplet K 2 . 



SEPARATED ATOM LIMIT 

A system of three well-separated atoms is more con- 
veniently described in terms of atom-atom distances ri 2 , 
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r23 and 7"3i and the angles 0i, <p 2 and 03 measured at 
each atom. The interaction energy for such a system is 
conventionally represented in terms of pairwise additive 
and nonadditive terms, 



V{v) =^2v diineI (r ij )+V 3 {r), 



(3) 



i<3 



where we use the shorthand (r) to indicate 
(f*i2, T"23, T3i, (pi, (p2, <t>?>)- For a pair of S-state atoms the 
long-range energy is 



Vdimor(r) 



-C 6 r~ 



C 8 r- 8 -C w r- 10 + 0(r- n ). (4) 



The nonadditive energy has several long-range contribu- 
tions that need careful consideration. The best-known 
is the Axilrod-Teller-Muto (ATM) triple-dipole (DDD) 
term, which has the form [24| 

F 3 DDD (r) = Z m ^iii(r), (5) 

where 

Win(r) = 3(1 + 3cos0i cos02cos0 3 )rf 2 3 r 2 ~ 3 3 r 3 ~ 1 3 . (6) 

The triple-dipole term is one of several that arise in third- 
order perturbation theory from terms of the type 



(000\H[ 2 \lmO)(lmO\H 23 \lOn)(lQn\H! sl \QOO) 
(AEf + AE? n )(AEl + AE*) : 



(7) 



where the ket \lmn) indicates a product wavefunction 
with atoms 1, 2 and 3 in states I, m and n respectively 



and AE\ is the excitation energy for state I of atom i. 
The interaction Hamiltonian H'^ is usually expanded at 
long range in terms of multipolc-multipole interactions. 
The triple-dipole term arises when each of the three op- 
erators HL is a dipole-dipole interaction of the form 



HUDD) = 



[l t .[L 3 ~ 3(/Xj.fy)(^.fy) 



(8) 



where fy is a unit vector along r^. The triple-dipole 
term will be referred to below as a third-order (3,3,3) 
contribution to indicate the powers of the three distances 
involved. 

There are additional third-order terms that arise when 
one or more of the dipole operators is replaced by a 
quadrupole or higher-order moment. The low-order 
terms and their resulting contributions for 3 identical 
atoms are 



DDQ 



(r) = Zii 2 Wii2(r) 



y 3 DQQ (r) = Z l22 W 122 {v) 
(r) = Zn 3 JVii3(r) 



DDO 



(9) 
(10) 

(11) 

^ 3 QQQ (r) = Z 222 W 222 (r), (12) 
where Q stands for quadrupole and O for octopole. The 
coefficients Zi t i 2 i 3 are related to polarizabilities of rank 
l\, l 2 and I3 at imaginary frequencies, 



Z lMa = (l/ir) a ll (iuj)a l2 (iLu)a t3 {iLu)(L;. (13) 
Jo 

The corresponding geometric factors are 



(15) 
(16) 



3 

W n2 (r) = — [9cos</!> 3 - 25cos3(fe + 6cos(0i - 2 )(3 + 5cos20 3 )]rf 2 3 r 23 4 r 31 4 + c.p. (14) 
15 

Wi 22 (r) = — [3(cos^i + 5 cos 5(?!>i) + 2Ocos(0 2 - ^)(1 - 3cos20i) 
o4 

+ 70 COS 2(02 — ^3) cos< ?H] r r2 4r 23 5r 31 4 + C 'P- 
5 

Wu 3 {r) = —[9 + 8cos20 3 - 49cos40 3 + 6cos(0i - 2 )(9cos0 3 + 7 cos S^r^r^r^ 5 + c.p. 
15 

W 7 222(r) — -—{—27 + 220 cos 0i cos 02 cos0 3 + 490cos2(/>i cos20 2 cos2</> 3 
128 

+ 175^082(0! -0 2 )+cos2(0 2 -0 3 )+cos2(0 3 -0 1 )]}r- 2 5 r 23 5 r- 5 , (17) 

where c.p. indicates summation of all cyclic permutations of indices. It may be noted that the multipole operator on 
each atom appears in two terms in Eq. J7J), so that changing one dipole operator into a higher-order one affects two 
of the inverse powers. 

In addition to these third-order terms, there are terms arising from fourth and higher-order perturbation theory that 
make important contributions to the long-range atom-diatom coefficients. In principle, any combination of interaction 
operators could produce a 4th-order term, although only those that involve excitations on all three atoms produce 
a contribution to the 3-body nonadditive energy. In addition, any odd-order operators (dipoles and octopoles) must 
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occur an even number of times on each atom to satisfy parity constraints. Thus there are terms such as 

(0Q0|ff{ 2 (PP)|lm0)(?m0|ff 23 (PP)|fap) (lnp\H^ 3 (DD)\lqO)(lqO\H{ 2 (DD)\000) 
(Ai?/ + AEDiAEf + AE 2 + AE^AEf + AEj) ' 

which gives rise to a (6,6,0) contribution, and 

(000|iJ( 2 (DD)|Zm0)(/m0|^ 3 (DD)|?np)(/np| J ff( 2 (QQ)|( ? 0p)((70p|iJ^(DD)|000) 
(A^ + AEl){AE} + AEl + AE^AE] + AE*) ' 

which gives rise to a (8,3,3) contribution. However, terms such as 

(000 1 H[ 2 (DP) \lrnO) (ImO | H' 23 (DP) | Inp) {Inp \ H[ 2 (DP ) qOp) (q0p\H^ (DP) |000) 
(A^J + AEl)(AEf + AEl + A^ 3 )(A£i + AEf) 



(18) 



(19) 



(20) 



are forbidden by parity, so that there is no (6,3,3) con- 
tribution. The fourth-order dipole energy ijl%|l has been 
given within a Prude oscillator model by Bade 

F 3 DDDD (r) = -gvV (1 + cos 2 fa) r^ + cp., (21) 

where V is a characteristic excitation energy and a is 
the static atomic dipole polarizability. If we adopt the 
definition Zuu = (5/32)V r a 4 , which corresponds more 
generally to 



ZhlilsU = (1/t) / ah ( iuj W 2 [iu))a h (ioj)a h (iu) do, 
Jo 

(22) 

then 



V, 



DDDD 



(r)=ZiuiWii U (r), 



(23) 



it within the Prude model of ref. |27]. We adopt the 
definition Z\xn\ = (35/256)T^a 5 , which corresponds to 



Ziim = (l/7r) / [a 1 ^)] 5 duj. 



where 



(25) 

In the 5th-order case, Eq. (10) of ref. |2j]then reduces to 
^DDDDD (r) = ZuiuWmuir), (26) 

with 

Winn (r) = 9 [2 + 3 cos 2 fa + 3 cos 2 fa + cos 2 fa 

+ 9 COS fa COS fa COS </>3] 7"i2 9r 23 3? "31 3 

+ c.p. (27) 



CONNECTING THE ATOM-DIATOM AND 
SEPARATED ATOM LIMITS 



Wiui (r) 



1 



cos 2 fa) r 12 6 r 13 6 



c.p. 



(24) 



In addition, there are fourth-order terms like Eq. 1181) 
but with one pair of the dipole-dipole operators replaced 
with dipole-quadrupole operators. These give terms with 
powers (8,6,0), (8,3,3), (7,7,0), (7,4,3) and (6,4,4). 
The angular factors have been given to within overall 
scaling factors by Lotrich and Szalewicz (though 
the (6, 3, 3) term that they describe is in fact zero and 
there are several typographical errors in their equations). 
There is also a fourth-order term involving a single 
octopole-dipole operator (with all the rest dipole-dipole) 
that also contributes to (8, 6, 0) but is not mentioned in 
ref.Hl 

There are also analogous terms arising from fifth-order 
perturbation theory, which are also constrained by the 
requirement that odd-order operators must occur an even 
number of times for each atom. The leading such term 
(PPPPP) has powers (9,3,3), but its explicit angular 
form has not been given previously. We have evaluated 



To make an explicit connection between the atom-atom 
additive and nonadditive dispersion terms and the atom- 
diatom dispersion formulae @ , we must collect terms in 
the (ri2, T23, r 3 i) representation that contribute to indi- 
vidual powers of the Jacobi distance R. Let us consider 
the case in which a diatom made up of atoms 1 and 2 
is separated from atom 3. For simplicity, we will assume 
that the three atoms are identical. The atom-atom dis- 
tances may be written in terms of Jacobi coordinates R, 
r and 9, 



r\2 



r23 



'";si 



R 2 + — + Rr cos ( 



R 2 + — - Rr cos ( 



(28) 
(29) 

(30) 



and the cosines of fa, fa and fa are given by the cosine 
rule. The approach we take is to express the different 
contributions to 3-body energies in Jacobi coordinates 
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using these equations and then to expand the results as 
power series in r/R. 

We consider first the additive terms Vdi mer ( r ij)- F° r 
R ^> r the atom-atom Cq term makes contributions to 
the atom-diatom coefficients of Eq. given by 



Cf dfi (r,9) = 2C 6 ; 



Cf dfi {r,9) = C 6 ( ~ + 8P 2 (cos9) ) r 2 ; I •!2 ) 



(31) 



4 7 



7 



(33) 



Similarly, the atom-atom Cs and Cio terms contribute 
C 8 add ' 8 M) = 2C 8 ; (34) 



C? o dd ' 8 (r,0) = C 8 f^ + ^P 2 (cos ( 



C 



add, 10 



(r,9) = 2C 



10- 



(35) 
(36) 



It may be seen that the atom-atom pair potential con- 
tributes no anisotropy to the longest-range (i? -6 ) term 
in the atom-diatom potential. All such anisotropy must 
come from 3-body nonadditive terms in the potential. 
The only third-order 3-body term that contributes to 
Cq(t,9) is the triple-dipole term. Its geometric factor 
may be expanded at large R, 



J 



Wiu(R,r,6) = -6P 2 (cos9)r- 3 R- 



6 



P 2 (cos( 



36 



P 4 (cos9) )r- 1 R- 8 + 0(rR~ 



10 \ 



2 7 v 7 

Similarly, the geometric factors for the third-order DDQ, DQQ and DDO terms may be expanded 

120 

Wn 2 {R,r,9) = —[P 2 (cos9)~P 4 (cos9)}r- 3 R- 8 + 0{r- 1 R- w ) 

W 122 (R,r,9) = 30P 4 (cos9)r- 5 R- 8 + 0(r- 3 R- 10 ) 
Wu 3 (R,r,9) = -40P 4 (cos9)r- 5 R- 8 + 0(r- 3 R- 10 ), 



(37) 



(38) 

(39) 
(40) 



where any contributions from cyclic permutations are now included. These three terms thus contribute to the atom- 
diatom Cs(r) and its anisotropy but not to Ce(r, 9). 

As noted above, the fourth-order DDDD term does not have a (6,3,3) contribution but does have a (6,6,0) contri- 
bution. It can thus contribute to the atom-diatom Ce coefficient. Its geometric factor (124(1 may be expanded at large 



Willi (-R, r,6) 



-6 r>-6 



6 [2 + P 2 (cos 9)] r~ b R 



33 402 

Y + - 



P 2 (cos( 



144 

T 



P 4 (cos<9) r- A R- & + 0{r- 2 R- 10 ). 



(41) 



Collecting these equations together provides expressions for the behaviour of the atom-diatom dispersion coefficients 
C„{r) at large r, 



C° 6 (r 
Clir 

C° s (r 
Cl(r 



2C 6 + 12Z nur - 



0(r~ 



&Z ln r- 3 + 6Z lllir - 6 + 0(r- s ); 

\c,r 2 + 2C 8 \Z^T~ X + f Z lm r- 4 + 0(r" 6 ); 

SC 6 r 2 + ~Z xxx r^ + ™Z 112 r- 3 + ^Z ml r- 4 + 0(r-% 



120 



144 



= — Z^r- 1 - —Z 112 r- 3 + —Z lllir ~ 4 + (40Z 113 - 30Z 122 ) r- 5 + 0(r~ e ). 



(42) 
(43) 

(44) 
(45) 
(46) 



The 0(r 8 ) terms in Eqs. f4*2*j) and lj4*3l come from the order energy, while the 0(r 6 ) terms in Eqs. to 
(8,6,0), (8,3,3) and (6,4,4) contributions to the fourth- come from both these and the (7,7,0) and (7,4,3) con- 
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tributions. The fifth-order term (Eq. does not con- 

tribute until 0(r~ 9 R~ 6 ). 



FITTING DISPERSION COEFFICIENTS FOR Li 
+ Li 2 

The equations above apply when all three of r-y, r 2 
and r 3 are large. However, when any of them is small, 
the power series is insufficient. We therefore use Eqs. 
(I42|l to H46|) as the long-range limits of more general ex- 
pressions, constructed by (i) multiplying the individual 
inverse power terms by damping functions D n (r) and (ii) 
adding a short-range exponential term to allow for the ef- 
fects of orbital overlap. 

As described above, Rerat and Bussery-Honvault |2l) 
have calculated isotropic and anisotropic dispersion co- 
efficients for Li + Li2 ( 3 S+) as a function of r and have 
fitted them to long-range expansions. They stated that 
"no asymptotic form of the coefficients C§ exists to our 
knowledge", but found empirically that Cf(r) required 
both r -3 and r~ 6 terms, while Cg (r) required only r -6 . 
Our expressions (|42|) and l|43[) above provide the expla- 
nation for this. However, our results also show that there 
should be relationships among the coefficients of the fit, 
and the resulting constraints were not included in ref. 
EH In particular, the coefficient of r~ 6 in Cf should be 



half that in Cg and Zim can be related at least approx- 
imately to Zin and C$ as described below. 

In devising functional forms for Cg(r) and Cf(r) it 
is important to consider damping of the inverse power 
terms. For the two-body interaction energy, the most 
popular approach is to use Tang-Toennies damping func- 
tions |30| of the form 

D n (R) = l-e-" R j2^r- ( 4? ) 

fe=0 

When damping 3-body terms such as Eqs. © and i|14|) 
to i|17|) . a damping function is required for each r™. We 
have chosen to use ^ D 2n (ri) rather than D n (ri) for this 
purpose, because this recovers the correct D n (R) in the 
2-body energies. Thus we damp 3-body terms according 
to prescriptions such as 

r^r^T-ii 3 -» r^r^rii 3 y/D 6 (r 12 )D 6 (r 23 )D e (r 31 ) ; (48) 
rgr^r^ - rgr£r£ V 'As (r 12 )D s (r 23 )D s (r 31 ) .(49) 

When damping is introduced in this way, the r -6 terms 
in the expressions for dispersion coefficients are damped 
with Dq(t), but the r -3 term in Cf (r) is damped with 
\J Dg(r) rather than D 3 {r). The expressions that we fit 
to are therefore 



C 6 °(r) = 2C 6 + 12Z lllir - 6 D 6 {r)+Aexp{-Cx); 



C 2 e (r) 



6Z llir ~ 3 y/D 6 (r) + 6Z lllir - 6 D 6 (r) + B exp(-Cx), 
I 



(50) 
(51) 



where x — (r — r )/r and r = 7.0 A. 

We have fitted the values of Cf(r) and Cf (r) in ref. 
EH to the func tional for ms lfS7)|) and JSU- [Our C*(r) 
are C^ 0A (r)/ v /A(A + 1) in the notation of ref. 21.] The 
triple-dipole coefficient was fixed at a value taken from 
variational calculations with Hylleraas basis sets by Yan 
et al. 31], Ziu = v ahc jZ = 5.687 x 10 4 E h al (where a 
is the Bohr radius and Eh is the Hartree energy). There 
is no ab initio value available for the quadruple-dipole 
coefficient, but in a Drude model Zim = (5/32)V^a 4 as 
described above, where V is a characteristic excitation 
energy and a is the atomic dipole polarizability. Within 
the same model, C 6 = (3/4)Va 2 and Z ln = (3/16)Va 3 . 
Combining these results gives an estimate 



Zi 



in 



lO(^m) 2 
3C fi 



(52) 



For Li3 with Cq and Zm values from ref. l3lL this gives 
Z 11U = 7.73 x 10 6 E h al 2 . The C 6 coefficient in Eq. ^5Djl 
could not be fixed at the value C§ = 1393 E^a^ from ref. 



Ell because the results for Ceir) in ref. EH converge on a 
slightly different asymptotic value, so C§ was allowed to 
vary in our fit. 

Meath and coworkers |U |3^| have calculated disper- 
sion damping functions for Li2 and found that the r~ e 
damping function is around 0.45 at the diatomic mini- 
mum, r e w 4.2 A. With the Tang-Toennies form of D e (r), 
this requires 6 k 1.5 A -1 . We found that the exponent 
C needed to fit the short-range part of C§(r) was very 
different from 6ro, so we allowed brg and C to be dif- 
ferent and fixed b at 1.5 A -1 . The remaining param- 
eters were then determined by a weighted least-squares 
fit, giving values C 6 = 1414.8 E h a 6 , A = 17.2 ±2.6 E h al, 
B = 35.2 ± 3.9 E h al and C = 5.13 ± 0.25. The quality 
of the resulting fit is shown in Fig. together with the 
contributions of the individual terms in Eqs. (|50|l and 

(50. 

Rerat and Bussery-Honvault [ilj have also given val- 
ues of Cg(r), Cf (r) and C|(r) for the single distance 
r = 4.2 A, which is near the equilibrium distance for 
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1 ~T-^-r-^ . =-^ 1 , ^—-^ ■ — 

2 4 6 8 10 12 14 

r/A 



FIG. 1: Fits to dispersion coefficients Cg(r) (upper panel) 
and Cf (r) (lower panel) for Li + Li2 interactions. 



triplet hi2- In our notation their values correspond to 
3.92 x 10 5 E h a%, 1.91 x 10 5 E h al and -0.38 x 10 5 E h a% 
respectively, where Eh is the Hartree energy and do is 
the Bohr radius. Evaluating the long-range expressions 
105 to gBJ} at r = 4.2 A, using C 6 = 1393 E h al and 
£m = 5.687 x 10 4 E h a% as above with C 8 = 83426 E h al 
H3, Zi 12 = 5.81 x 10 5 E h al\Zi22 = 6.41 x 10 6 E h al 3 
and Z U3 = 1.70 x 10 7 E h a\ 3 |3, gives 4.076 x 10 5 E h a% 
8.402 x 10 5 E h al and +0.725 x 10 5 E h al, respectively. 
It thus appears that r = 4.2 A is too small a distance 
for the Cs (r) coefficients to be dominated by undamped 
long-range contributions, and indeed it may be seen from 
Fig- El that the exponential contributions to Cg(r) and 
C|(r) are quite large at this distance. It would be very 
interesting to calculate r-dependent Cs coefficients and 
explore the onset of long-range behaviour. 



GLOBAL FUNCTIONAL FORMS FOR 
POTENTIAL ENERGY SURFACES 

The results obtained above have important implica- 
tions for the choice of functional forms to represent 
potential energy surfaces for triatomic systems. For 
low-energy scattering calculations, it is highly desirable 
to have an interaction potential that dissociates cor- 
rectly both to three atoms and to an atom and a di- 
atomic molecule. This is especially important for pro- 
cesses such as 3-body recombination and collisions of 
long-range diatomic molecules with atoms, which are of 
current interest in studies of cold molecule formation 

i i a 1 i m mm in la e in m 

A global potential energy surface for a homonuclear 
triatomic system must be symmetric in the atom indices 
if it is to reproduce the full symmetry of the system. 
The simplest coordinate system that achieves this is a 
set of 3 bond length coordinates (^12, ^23, ^3i)- Jacobi 
coordinates cannot easily describe the full symmetry. 

Our approach is to obtain a potential energy surface 
from high-quality ab initio calculations on a grid of val- 
ues (^12, r23, ?"3i). Once this has been done, we need to 
interpolate and extrapolate it in a way that incorporates 
the correct long-range behaviour. In one dimension, re- 
ciprocal power - reprodu cing kernel Hilbert space (RP- 
RKHS) interpolation j^l [3I [37J provides an attractive 
way to obtain a potential with desired inverse power be- 
haviour at long ran ge. However, multidimensional RP- 
RKHS interpolation |l8j,|35| at best gives a potential that 
extrapolates beyond the points as a simple product of 
inverse powers in the different coordinates. None of the 
long-range 3-body terms have this form, so a more so- 
phisticated approach is required. 

Two different cases may be distinguished. For systems 
such as spin-polarized Na3 or K3, the nonadditive forces 
are substantial at short range but do not dwarf the addi- 
tive forces 0,E|. Under these circumstances, we have 
found it best to transform the potential to a form that 
does behave as a simple product of inverse powers at 
long range and then interpolate in that form. The first 
stage in this process is to subtract the (assumed known) 
sum of pair potentials from the ab initio points accord- 
ing to Eq. (0 to obtain the nonadditive contribution to 
the interaction, V^r). The leading terms in the long- 
range expansion of Vs(r) arc the DDD and DDQ terms. 
However, it may be noted that the DDD term of Eq. © 
vanishes on a seam in the angular space and the DDQ 
term of Eq. I|14l) vanishes at all linear configurations. It 
is therefore not adequate to divide V3(r) by an angular 
factor in order to obtain a form that depends only on in- 
verse powers. Fortunately, the coefficients Z\\\ and Z\\2 
are known for all the homonuclear alkali metal atom sys- 
tems [34| , so that damped versions of these terms can be 
subtracted from the total nonadditive energy V3 to give 
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a quantity V 3 , 

V s '(r) = Vs(r) 



^3 D d°mp( r ) 



^3, damp \ l ) 



(53) 



In our work on K 3 |l5| , we used a single damping function 
for both terms. However, the present work has shown 
that it would be better to choose separate damping func- 
tions for each inverse power term as in Eqs. (|48|l and (|49ll . 
and to define (for example) 



DDQ 



3, dam 



p (r) = -tZii 2 [9cos( 



, 25cos3c/> 3 

16 

6cos(</)i — <f> 2 ) (3 + 5cos2^ 3 )] 



V ' D&{ri 2 )Dg,{r 23 )D & {r 3l ) 

X n 4 A T C.p. 



r 12 r 23 r 31 



(54) 



The leading term in the asymptotic multipole expan- 
sion of V 3 is the fourth-order dipole-dipole-dipole term 
(DDDD), which has the more complicated (unfactoriz- 
able) form of Eq. (f2~4"j) above, with powers (6,6,0). If Zim 
is known, this too could be subtracted out. However, this 
term is negative at all geometries, so a satisfactory alter- 
native is to eliminate it by defining Vg"(r) = g(r) x ^'(r), 
where 



12' 23' 13 



(1 + cos 2 <j)i) r^ 3 + (1 + cos 2 <j) 2 ) rf 3 + (1 + cos 2 cj> 3 ) r\ 



(55) 
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The leading asymptotic term of the function V 3 " now has 
the form —constant x r^ 2 3 r 23 ri 3 and is suitable for an 
"isotropic" extrapolation of the type that results from a 
multidimensional RP-RKHS interpolation. 

The approach that we use is therefore to construct V 3 " 
at the ab initio points as above and then interpolate it 
using the fully symmetrized 3D RP-RKHS interpolation 
method [l^. An RP-RKHS interpolation with respect 
to r p with RKHS parameters n and m gives a potential 
with leading long-range powers r~ p ( m+1 ) and r -p( m + 2 ). 
Appropriate choices thus include p — 3, m — (as in our 
work on K3 [l5j) and p = l,m = 2. The interpolated 
potential is then rebuilt as 



Vs(r) 



S(r) 



T^DDD ( \ , y 



DDQ 
3, damp 



(r) 



(56) 



A system such as quartet Li3 requires a different ap- 
proach 01- I n this case the nonadditive forces are so 
large that at short range it does not make sense to decom- 
pose the potential into additive and nonadditive parts at 
all. The nonadditive potential is many times larger than 
the additive potential [lflj l , and the decomposition would 
require the final short-range potential to be expressed 
as the difference of two large numbers. Nevertheless, at 
long range a decomposition according to Eq. is es- 
sential. Under these circumstances, we found it best to 
carry out an unconstrained fit to the ab initio points at 
short range, without imposing the correct long-range be- 
haviour, to obtain a short-range function Vsr(i"). Wc 
then use a switching function S(r) to join this onto the 
correct long-range form. We thus have 

V(r) = S(r)V SK (r) + [1 - S(v)]V LR (r). (57) 

The long-range form must be valid when any of the atom- 



atom distances is large. In our work on Li3 [itI ]. we used 
V w {r) = V V dimel ( rij ) + U 3 , LR (r), (58) 



where 



V3,LR.(r) = ^p( r )+CaSp( r )+^?(r)+V^? P (r). 

(59) 

The function of the term V^^^r) is to ensure that 
the atom-diatom dispersion coefficients Cg(r) and C|(r) 
have the correct values (given by Eqs. (|5T7|l and (|5*T|l % ) even 
when one of the atom-atom distances is small. When 
R^> r, this is achieved by defining 

CW - -\A + BP 2 (cos 0)} 



x e^p{-Cx)r 23 3 r 31 3 y/D 6 (r 2 3)D 6 {r 31 ) 
+ c.p., (60) 

where x = {r\ 2 — ro)/ro as before and the parameters A, 
B and C come from fits to numerical values of Cg (r) and 
C|(r) as described above. In evaluating Eq. H60|l it is 
convenient to use an approximate form of P 2 (cos 9) that 
is valid for R ^> r but is well-behaved at all geometries, 



P 2 (cos( 



1 



(1+3 cos 4>i cos 4» 2 cos fa) 



(61) 



This is already evaluated as part of Wm(r), and since 
it is symmetric the cyclic permutations required in Eq. 
(|60|l involve no extra geometric calculations. 

Finally, the switching function S(r) must become zero 
when any of the three atom-atom distances is large. For 
Li 3 we chose to use 



1 



S(r) = - tanh[l - si(?'i + r 2 + r 3 - s 2 ) 



(62) 
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with the parameters s\ and S2 determined in such a way 
that the switching takes place in a region where both 
functional forms give reasonably accurate energies. 



CONCLUSIONS 

We have investigated the relationship between long- 
range intermolecular forces for triatomic systems in the 
atom-diatom and atom-atom-atom representations. We 
have obtained expressions relating the dispersion coeffi- 
cients in the two representations. We have shown that 
the anisotropy of the atom-diatom Cq dispersion coeffi- 
cient arises entirely from nonadditive terms in the 3-body 
expansion. The most significant contributions at long 
range arise from the third-order triple-dipole term and 
the fourth-order quadruple-dipole term. The leading con- 
tributions to the atom-diatom C% coefficient arise from 
the additive atom-atom Cg and Cs coefficients and the 
third-order nonadditive triple-dipole and dipole-dipole- 
quadrupole coefficients. 

There is great current interest in the formation of di- 
atomic molecules in cold atomic gases, and the collisional 
properties of such molecules are of great importance. 
Calculations on these collisions need triatomic interac- 
tion potentials that dissociate properly both to an atom 
and a diatomic molecule and to three separated atoms. 
We have used our results to suggest strategies for obtain- 
ing such potentials. 
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